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Abstract 

The analytical calculation of the self-referred lacunarity is used as a validation standard of the 



computational algorithm. In this supplementary material to our article (see cond-mat / 0407079( ) 



we present a detailed calculation for two simple shapes, namely a square box and a cross. 

1 Introduction: self-referred lacunarity 

The lacunarity of an object is defined as 

where Z^^^ and Z^^^ are the first and second moments of a mass distribution and (T^(r) and /x(r) 
are its variance and mean. We can calculate fi{r) and (T^(r) by defining a function A{x,y,r) that 
gives the area of the object inside a window of radius r, whose center is in the point in the 

following way 

IxIyMx,y,r)dxdy 



Ix ly dxdy 



^ (r) = — rr-TTi • 2 

L ly dxdy 

For the self-referred approach of lacunarity the integral is calculated only for values of x and y that 
fall over the object, see Figure ^ 

2 Analytical calculation for simple shapes 
2.1 Object: a square box 

We first show how to determine the self-referred lacunarity for a solid square of side L shown in 



Figure 1(a) In order to calculate the analytical lacunarity of a square it is necessary to determine 



the area funcition A{x,y,r). This function is determined by dividing the square in regions labeled 



by /, II and /// as in Figure 1(a) and 1(b) There is one region of type / and four of type // and 



///. It is also necessary to divide this calculation into two intervals: r = [0,L/2] and r = [L/2,L]. 
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Figure 1: Square regions that correspond to the self-referred lacunarity calculation for the intervals (a) r G [0, L/2\ 
and (b) r £ [L/2,L\. The sliding windows in the region I are completely filled in (a) and capture all the 
square in (b). In the regions // and /// they are partially filled in (a) as well as in (b). 



For the first interval, r € [0, L/2], the functions A\{x,y,r) and respective integration ranges for 
X and y are 

y, r) = 4r^, x € [r, (L — r)], G [r, (L — r)] 
A]j{x,y,r) = 2r^ + 2ry, x G [r, (L - r)], G[0,r] 
Alli^^y^r) = r"^ + [r + {L-y)]x + r{L-y), 
x€[0,r], ?/G[(L-r),r] 

The mean value is determined by the integration of the A\ over their respective regions divided by 
the total area of the square, according the Equation ^ 

!xIyT.i Al{x,y,r)dxdy 



Mr) 



L ly dxdy 



J2 J J A]ix^y^r) + 4:A]j{x,y,r) + 4A]jj{x,y,r)dxdy 



L2 

The four regions II and /// are equivalent and it is necessary only to integrate over one of region 
II and /// and multiply the result by 4. The variance is calculated in analogous way according to 
Equation |21 



= J{A]{x,y,r)-f,^{r))' 
+ A{A]j{x,y,r) - fii{r))^ 

+ 4.{A]jj{x, y, r) - fii{r)fdxdy 
_ r^{r -6L){3r - AL){3r -2L) 
9IA 



(3) 



2 



interval 



[/ii(r)]2 9(r-2L)4 

The calculation for the second interval, r € L], proceeds in the same way, with the Af functions 
and their respectives integration intervals given by 



A'j{x,y,r) 



X G [(L — r), r], y £ [{L — r), r] 



Ajj{x,y,r) = L{r + y), xe[{L-r),r], y G [0, (L - r) 
A'jii{x,y,r) = r"^ + [r + {L-y)]x + r{L-y), 
xG[0,(L-r)], ye[r,L] 



and the mean, variance and Lacunarity as in 

r2(r - 2L)2 



A2(r) 



L2 

(3L - r)(r - L)^{3r^ - Ur^L + 12r2L2 + 6rL^ - 



9L4 



3\2 



9r4(r - 2L)4 

The final expression for the whole interval self-referred lacunarity is 

r G [0,L/2]; 



4L^(5r-6L)^ 
9(r-2L)4 ' 



A(r) 



L\2r^-GrL^+L^f r G /2 M 



Figure 12 shows a plot of the expression above for a square of length L = 100, in absolute units of 
length. 
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Figure 2: The analytically calculated self-referred lacunarity curve for a square with side L of 100 absolute units of 
length. 
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The analytical procedure to calculate the lacunarity of a cross, which occupies a region k, is similar 
to the procedure presented above for a square. Again we define a function A{x,y,r). For this case, 
there are seven intervals and each one was divided into regions which define an associated function 
Ai{x,y,r). Figure ISl shows the first interval considered in the self-referred lacunarity calculation of 
a cross with dimensions LI and L2 and the respective regions. Figure 01 presents only non-repeated 
regions i, however the calculation regards all the possible regions j. In the following we show the 
calculation of the first interval in more detail. For the remaining intervals we provide only the 
functions Ai{x,y,r), (T^{r) and A(r). Figures 0] to [HI presents the divided regions for the 

remaining intervals labeled from 2 to 7. The expression for lacunarity, A(r) = (cj^(r)//i^(r)) + 1 
where ^ and cr^ is given by 



H{r) 



Ikdk 



and 



Ik [EAj{x,y,r) - n{r)fdk 
Ikdk 

where dk is the total k area {Ak = 4L1L2 + L2^) of the cross. 



2.2.1 Interval < r < ^ 
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Figure 3: The cross divided to regions for the first intervaL 



Ai{x,y,r) = xe[r,L2-r], y e [-{L2 - 2r),0] 

Aii{x,y,r) = 4r^ xe[r,L2-r], y G [0, LI] 
Ajii{x,y,r) = {2r + Ll-y)2r xG[r,L2-r], yG [LI, Ll+r] 
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Avix,y,r) = (r + x)(2r) xG[0,r], y G [2r, LI] 
Avi{x,y,r) = Ar^ — y{r — x) xS[0,r], y G [r, 2r] 
Avil{x,y,r) = Ar'^-y{r-x) xG[0,r], y G [0, r] 

The Calculation of the mean /x for the first region gives 

+ AAii + AAiii + 8yljy + 8Av + ^Ayi + 4Avii)dk 



Ak 

r^{4L2^ - 4rL2 + Sr^ + 16L1L2 - 8rLl) 



L2(4L1 + L2) 

We make LI = 2L e L2 = L to simplify the expressions. Thus we have for the previous expression 

, , r^(36f -20Lr + 3r^) 
f^dr) = • 

In an analogous way cj^(r) is given by 

([Ai - ^1I{L1, L2, r)f + ... + [iAvii - /i/(Ll, L2, r)f) dk 



aj{r) 



Ak 

384L2L12 + 576rLl2 - 432r2Ll - 288L1L22 



9L22(4L1 + L2)2 
+496rLlL2 - 216L2r2 - 48L2^ + 81r^ + 124rL22)r^ 

Simplifying the above expression we have 

2 (9r3 - 120Lr2 + 380rL2 - 240L3) 

^^^"^ - sTl^ • 

The lacunarity expression in this interval is finally given by 

. , , 4L2(-300Lr + 59r2 + 324L2) 
A/(r) = 



(36L2 - 20Lr + 3r2)2 



2.2.2 Interval ^ < r < L2 

Ai{x,y,r) = 2{2rL2)-L2^ x G [L2 - r, r], ?/ G [-(L2 + 2r), 0] 
An{x,y,r) = {2r + y - L2)L2 + {L2 - y)2r x£[L2-r,r], y G [0, L2] 
Ajii{x,y,r) = L2x2r xe[L2-r,r], y G [L2, LI + L2 - 2r] 
Ajv{x,y,r) = L2{Ll + L2-y) xe[L2-r,r], 

y e[Ll + L2- 2r, LI + L2 - r] 
Av{x,y,r) = {LI + L2 - y){r + x) xG[0, L2-r], 

y e[Ll + L2- 2r, LI + L2 - r] 
Avi{x,y,r) = 2r{r + x) xG[0,L2-r], y G [L2, LI + L2 - 2r] 
Avil{x,y,r) = 2r{r + x) + {r - x){L2 - y) x G [0, L2 - r], y G [L2 - r, L2] 
Avin{x,y,r) = 4r^ - (2r + y - L2)(r - x) xG[0,L2-r], yG[0,L2-r] 
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Figure 4: The cross regions to interval ^ < r < L2. 



r2(-20rL + 36L2 + 3r2^ 
91? 



1 



-{l2f)r'L - UL'r - 9r« + 1248r^L^ + 452LV^ 



- - GOAr^L^ - 1128LV^) 



A//(r) 



(-12 L^r - 192 r^L + 452 L^r^ + - 36 L^r^ + 12 + 168 r^'L^) 
r4(-20rL + 36L2 + 3r2)2 



2.2.3 Interval L2 < r < ^^^^ 



Aj{x,y,r 
Aji{x,y,r 
Ajii{x,y,r 
Aiv{x,y,r 
Av{x,y,r 



2L2(2r) -12"^ x G [0, L2], ye [-L2, 0] 

4rL2 -L2^ x e [0, L2], ye [0, r - L2] 

L2{r + y) + {r - y)2r xe[0,L2], y e [r - L2, LI - r] 

{LIL2) + {r - y)2r xG[0,L2], ye[Ll-r,r)] 

{Ll-{y-r))L2 xG[0,L2], y G [r, LI] 



, , 4L ^ In 2 2 
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Figure 5: The cross regions for the interval L2 < r < LLl±i^^ 



^^^^ ^ 27 9 9 27 81 



. , , 64 + 320 r^L^ - 348 Lr^ - 204 L^r + 37 
Am(r) = -3 



(12rL + L2 + 6r2)2 



2.2.4 Interval < r < LI 



Ai{x,y,r 
Aji{x,y,r 
Ajii{x,y,r 
Aiv{x,y,r 
Av{x,y,r 



2L2(2r)-L22 xG[0,L2], ye[-L2,0] 

2rL2+{r-y-L2)L2+{y + r)L2 xe[0,L2], ye[0,Ll-r] 
L1L2 + 2rL2 + {r - y - L2)L2 x G [0, L2], y G [LI -r,r- L2] 
LlL2+{r-y)2r xG[0,L2], y G [r - L2, r)] 
(Ll-(y-r))L2 xG[0,L2], y G [r, LI] 



/ N 4L ^ 1 2 2 2 



4 4 68 2.2 20 , , 236 212 . 

9 9 27 27 81 



A/y(r) 



L (256 r^L + 28 - 228 rL^ + 71 L^) 
(12rL + L2 + 6r2)2 
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Figure 6: The cross regions for the interval (-^^+-^^) < r < LI. 



2.2.5 Interval LI < r < LI + ^ 
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Figure 7: The cross regions for the interval LI < r < LI + —. 
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Ai{x,y,r) = 2(2rL2)-L2^ x € [r - LI, LI + L2 - r], 

y G [-L1 - L2 + r,-r + LI] 
Aji{x,y,r) = LlL2 + 2rL2 + {r - y - L2)L2 x e [r - LI, L2 + LI - r], 

y G [-r + Ll,r - L2] 
Ajii{x,y,r) = LIL2) + {r - y)2r x £ [r - LI, L2 + LI - r], ye[r-L2,Ll] 
Aiv{x,y,r) = L1L2+ {r - y){r + X + LI) a;G[0,r-Ll], y G [r - L2, LI] 
Av{x,y,r) = LlL2 + {r + x + Ll)L2 + {r-y-L2)L2 x£[0,r-Ll], 
ye[0,r- L2] 

Avi{x,y,r) = LlL2+{r + x + Ll)L2 + {r-y-L2)L2 x£[0,r-Ll], 
ye[-{r-Ll),0] 

, , -100 L^r + 90 r^L^ - 24 Lr^ + 2 + 49 
Mr) = 



Ay(r) 



a^{r) = - (398 + 6704 LV2 + 5717 L^r^ + 448 r^^L^ 
-2144 r^L^ - 8516 L^r^ - 2542 L^r - 48 r'^L + 2 r^) 

L2 (1605 L^ + 5412 L^r^ + 1662 r^L^ + 40 - 432 r^L - 3320 L^rS _ 4716 L^r) 
(-100 + 90 r2L2 - 24 Lr^ + 2 + 49 L^)^ 



2.2.6 Interval LI + ^ < r < LI + L2 

Ai{x,y,r) = 4L1L2 + L22 x G [LI + L2 - r, r - LI], 

y G [-(r-Ll),-(Ll + L2-r)] 
Aii{x,y,r) = L2^ + 3LIL2 + {r -y - L2)L2 x e [-{L2 + LI - r),r - LI], 

y G [-(Ll + L2-r),r-L2] 
Ainix,y,r) = L1L2 + (r - y)(2Ll + L2) x G [L2 + LI - r, r - LI], 
y G [r-L2,Ll] 

Ajv{x,y,r) = L1L2 + {LI + r + x){r - y) x G [0, L2 + LI - r], y G [r - L2, LI] 
Av{x,y,r) = LIL2 + (r - y - L2)L2 + {LI + r + x)L2 x G [0, L2 + LI - r], 
y G [0, r - L2] 

Avi{x,y,r) = LIL2 + {r - y - L2)L2 + {LI + r + x)L2 x G [0, L2 + LI - r], 
y G [-(Ll + L2-r),0] 



9 



LI 

(r-L2) 



IV 







VI 



-(Ll+L2-r) 

-(r-Ll) 



III 



Figure 8: The cross regions for the interval LI + ^ < r < LI + L2. 



Hvi{r) 



-100 L^r + 90 r^L^ - 24 Lr^ + 2 r"^ + 49 
9l2 



(11804 LV^ + 6881 + 464 r^L^ - 2360 r^L^ 



81L4 



-11766 - 7M2 L^r - 48 r^L + 2r^ + 2273 L^) 

L2 (4788 L^r^ + 666 r^L^ - 8 - 3180 L^r^ - 4284 + 2145 L'^ 
(-100 L^r + 90 r2L2 _ 24 Lr^ + 2 + 49 L'^f 



2.2.7 Interval LI + L2 < r < 2L1 + L2 

Ai{x,y,r) = 4L1L2 + L22 xG[0,L2], ye[-L2,0] 
An{x,y,r) = ALIL2 + L2'^ 2;G[0,L2], y G [0, r - LI - L2] 
Ani{x,y,r) = 3L1L2 + L2^ + {r - y - L2)L2 xe[0,L2], y G [r - LI - L2, LI] 



, , 31 2 20 ^ ^2 
fJ'Vll[r) = —L + — rL - -r 



2 
9 



A{2L-r){bL-rY 
81 
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Figure 9: The cross regions for the interval LI + L2 < r < 2L1 + L2. 



^Vllir) 



3L (13 - 780 rL^ + 4 + 48 r^L) 
(31L2 + 20rL-2r2)2 



Sumarizing the results for the cross, the whole expression for the lacunarity of a cross is given 
the following expression 



A(r) 



Aj{r] 

^11 {r) = 
Kin{r) -- 

Kiv{r) = 

Ay(r) = 

^vi{r) = 

^Vll{r) 



4L2(-300LT-+59r2+324L2) 
(36L^-20Lr+3r2)2 



2(-12L5r-192 r^L+452 r'^ +12 +168 r'^L'^ 

r4(-20rL+36 L'^+Sr'^)'^ 

_o 64 r^+320 r^L^ -348 Lr^-204 L3r+37 L^^ 
(12rL+L2+6r2)^ 

L(256r2L+28r3-228rL2+71L3) 
(12rL+L2+6r2)2 



(1605 L''+5412 L^r^+ie62 r'^L^+iO r''-432 r5L-3320 L3r3-4716 L^r) 
(-100 L^T+90 r2L2-24 Lr3+2 r4+49 L*)^ 

L2(4788 L''r2+666 r^L^ -8 r'^ -3180 L^r^ -4284 L^r+2145 L^) 
" (-100 L3r+90 r2L2-24 Lr3+2 r*+49 L^f 



3L(l3L^-780rL2+4r3+48r2L) 



(31L2+20rL-2r2)^ 

The corresponding intervals, I- VII, are shown in Tahel lTTTl and the resulting lacunarity curves 
this piecewisely determined function is plotted in Figure 
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Region 


Interval 


I 


r r O n 

re [0,f ] 


II 


r r O -1 

r £ [^,L2] 


III 


r € [L2, ^1+^2] 


IV 


r G [i^^i^,Ll] 


V 


r £ [LI, LI + ^] 


VI 


r G [LI + ^,L1 + L2] 


VII 


re [LI + L2,2L1 + L2] 



Table 1: Regions for the analytical calculation of the self-referred lacunarity. 




Figure 10: The analytically calculated self-referred lacunarity curve for a cross with dimensions LI and L2 of 100 
and 50 units of length respectively. 
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